Transport properties of directed percolation clusters at the upper critical dimension 
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We study the transport properties of directed percolation clusters at the upper critical dimension 
dc = 4 + 1, where critical fluctuations induce logarithmic corrections to the leading (mean-field) 
scaling behavior. Employing field theory and renormalization group methods we calculate these 
logarithmic corrections up to and including the next to leading correction for a variety of observables, 
viz. the connectivity, i.e., the probability that two given points are connected, the average two- 
point resistance and some of the fractal masses describing percolation clusters. Furthermore, we 
study logarithmic corrections for the multifractal moments of the current distribution on directed 
percolation clusters. 
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I. INTRODUCTION 

Directed percolation (DP) 0,0 is perhaps the simplest 
model for directed connectivity in disordered systems. It 
differs from conventional isotropic percolation (IP) 3] in 
that activity (in the following electric current) can perco- 
late only along a certain distinguished direction. A sim- 
ple and intuitive realization of DP is the random diode 
network (RDN), where nearest neighboring bonds on a 
tilted hypercubic lattice are randomly occupied (with 
probability p) with diodes, see Fig. ^ Evidently, the 
effective conductivity along the distinguished direction 
depends on the occupation probability p. If p is small, 
all cluster of connected sites are finite so that the con- 
ductivity vanishes in the limit of large distances. On 
the other hand, if p exceeds a certain critical threshold 
Pc, there is a finite probability to find a directed path 
between two points, leading to a finite resistance when 
averaged over many independent samples. At the criti- 
cal threshold p — Pc, the system undergoes a continuous 
phase transition where clusters of conducting paths dis- 
play a fractal structure. In contrast to IP, these clusters 
are anisotropic and at criticality they are rather self-affine 
than self-similar. 

In this paper we are interested in the electrical trans- 
port properties of DP clusters as realized by the RDN. 
Due to the self-affinity of the critical clusters, transport 
quantities like the average resistance M/{(x, x') between 
two connected points x and x' exhibit a scaling behav- 
ior as it is typical for (anisotropic) critical phenomena. 
Below the upper critical dimension dj, = 4 -I- 1 this scal- 
ing behavior is characterized by critical exponents which 
are, due to the effect of fiuctuations, anomalous. Above 
dc fluctuations are unimportant and the scaling behavior 
is purely of mean-field type. Right at dc, fluctuations 
lead to logarithmic corrections to the mean-field behav- 
ior. To date, the critical exponents for variety of trans- 
port quantities characterizing DP clusters are known to 
second order in an e expansion about dc, namely the 
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FIG. 1: (H-l)-dimensional directed bond percolation real- 
ized as a ransom diode network on a tilted square lattice. 
Electrical current can percolate only along the distinguished 
direction t. 



exponents for Mfl(x, x') 0, 
AfB(x,x'), Mred(x,x') and 



for the fractal masses 
n(x, x') respectively of 



the backbone, the red (singlyconnected) bonds and the 
chemical (shortest) path [j, Q as well as the exponents 
for the multifractal moments Mj'''(x, x') of the current 
distribution 0,0 (cf. Fig. Logarithmic corrections 
to these quantities have not been studied hitherto. With 
the computer resources and elaborate algorithms avail- 
able today, computer simulations can explore high di- 
mensions and it seems to be within reach to investigate 
the transport properties of DP clusters at dc numerically 
with good precision. Thus, we feel that it becomes im- 
portant to investigate the corresponding logarithmic cor- 
rections analytically. 

The system in statistical physics, for which logarithmic 
corrections have been studied perhaps most thoroughly, 
are linear polymers. The work of Schafer and cowork- 
ers P, showed that knowing the leading logarithmic 
corrections is, at least in the case of linear polymers, 
not sufficient to obtain a satisfactory agreement between 
theory and simulations. Rather, one has to push the 
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FIG. 2: Two-port setup for measuring the resistance of a DP 
cluster. The figure shows a generic current-carrying backbone 
when a current / is inserted at a point x and withdrawn at a 
point x'. The thickness of the bonds indicates the intensity of 
the current. The thickest bonds are the so-called red bonds 
that carry the full current. 



theory beyond the leading corrections, which involves, in 
addition to the basic analysis of the RG mapping, the 
calculation of scaling functions. We expect that higher- 
order logarithmic corrections will be likewise important 
for the transport properties of DP clusters. In part, this 
expectation is corroborated by the experience one has 
with logarithmic corrections in other percolation prob- 
lems. A comparison between numerical '11^ and ana- 
lytical ^3 results for dynamic isotropic percolation at 
the respective upper critical dimension clearly indicates 
the importance of the higher order logarithmic correc- 
tions. Simulations [13ij on dynamical aspects of DP at dc, 
on the other hand, show a satisfactory agreement with 
our theoretic predictions 0| already at the level of the 
leading logarithmic correction. There are no simulations 
available yet for a comparison to our analytic results on 
logarithmic corrections for transport on IP clusters p^ . 

In this paper we investigate the logarithmic corrections 
for transport on DP clusters up to and including the next 
to leading logarithmic correction. For simplicity, we as- 
sume that X = (0,0) and x' = (0,t). We calculate the 
quantities of interest as functions of the longitudinal or 
time- like variable t, that is Mii{t), Alsit), Mrod(i) and 
Mmin(i) as well as Mj^\t). By and large as a byproduct, 
we also obtain the connectivity, i.e., the probability that 
the points x — (0, 0) and x' = (0, t) are connected, P{t). 

The outline of this paper is as follows. In Sec. ^ 
we briefly review a generalization of the RDN, the ran- 
dom resistor diode network (RRDN), that certain ad- 
vantages over the RDN with respect to setting up a field 
theoretic model. Then, we explain our field theoretic 
model, including its variants, and highlight its physical 
contents. Section IIIII reviews rcnormalization group re- 
sults for the RRDN in order to establish notation and to 
provide known results that we need as an input as we go 
along. Section Hvl represents the main part of this paper. 
Here we derive the desired logarithmic corrections and 



we state our final results. Concluding remarks are given 
in Sec. [V] Technical details on the calculation of scaling 
functions are relegated to an Appendix. 

II. THE MODEL, ITS VARIANTS AND THE 
PHYSICAL CONTENTS 

A. The random resistor diode network 

The RRDN introduced by Redner [H [H [H is a 
simple model for electric transport in irregular media 
that comprises both, DP and IP. A RRDN consists of 
a c?-dimensional hyper-cubic lattice in which nearest- 
neighbor sites are connected by a resistor, a positive 
diode (conducting only in a preferred direction), a neg- 
ative diode (conducting only opposite to the preferred 
direction), or an insulator with respective probabilities 
p, P-, and q = 1 — p — p+ — p_. To be specific, 
let us choose n = l/Vd{l, . . . , 1) as the preferred direc- 
tion and assume that the bonds between two near- 
est neighboring lattice sites i and j are directed so that 
• n > 0. Moreover, let us suppose that the directed 
bonds obey the generalized Ohm's law 

-^..) (^.. J ^<..> - 4,,, ' (2.1) 

where VJ,^. = Vj — Vi is the voltage drop over the bond 
between sites j and i and Ih^, ,^ is the associated current. 
The bond conductances Ub = 0-76 are random variables 
where 7& takes on the values 1, {V), 9 {—V) (as usual, 
9 denotes the Heaviside function), and with respective 
probabilities p, p+, p_, and q, and where cr is a posi- 
tive constant. Note that this assumption means that the 
diodes are idealized, i.e., under forward-bias voltage they 
behave as "ohmic" resistors whereas they are insulating 
under backward-bias voltage. Further below we will also 
deal with nonlinear voltage-current characteristics under 
forward-bias of the type V ^ I"^ . The three dimensional 
phase diagram, i.e., the tetrahedron spanned by the four 
probabilities, features a non-percolating and three per- 
colating phases, viz. isotropic, positively directed, and 
negatively directed, with continuous transitions between 
the four phases. For a detailed discussion of the phase 
diagram see Refs. [iGlllTj . 

The long- length-scale behavior of the critical RDN and 
that of the RRDN at either of the 2 transitions from 
the non-percolating to the directed percolating phases 
are equivalent, provided of course, that the distinguished 
direction of the RDN corresponds to the distinguished 
direction of the RRDN or its opposite direction, respec- 
tively. Thus, we may investigate the transport properties 
of DP clusters as realized by the RDN via studying the 
RRDN at either of the transition the 2 transitions from 
the non-percolating to the directed percolating phases. 
We prefer to work with the RRDN because it has certain 
advantages with respect to setting up a field theoretic 
model. For details on deriving a field theoretic model 
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for the RRDN we refer to Refs. In the vicinity of 

the transitions from the non-percolating to either of the 
directed percolating phases this model can be written in 
the form of a dynamic response functional 0, 120I l2ll| , 



+ [0(Ao)-0(-Ao)] 



dt 
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A,A',A+A'/0 



(2.2) 



whose ingredients have a meaning as follows, — 
s^(x) is an order parameter field which lives on the d- 
dimensional real space with coordinates x = (r, t), where 
t ^ = X • n is the longitudinal part of x along the 
preferred direction and r = x_l is the corresponding 
((i_L — d — l)-dimensional transversal part. In addition to 
the space coordinate, the order parameter field depends 
on a Z?-fold replicated current variable A living on a D- 
dimensional replica space which we will explain a little 
further below. Regarding the current variable the order 
parameter field satisfies the constraint Sj5(x) = 0. The 
parameter r specifies the distance in phase space to the 
transition from the non-percolating to the directed phase 
of interest which occurs in mean field theory at r = 0. 
w is proportional to and p is a kinetic coefficient. 
The meaning of Aq will become clear shortly when we 
elaborate on the replica space. 

For regularization purposes, the replica space resem- 
bles a discretized D-dimensional torus, i. e., A = fcAA 
where fc is a D-dimensional vector with integer compo- 
nents fc(") satisfying -M < fc(") < M and k^"'^ = fc(") 
mod (2M). To extract physical quantities from the 
replica formulation one has to study the limit D — > 0, 
M 00 with (2A/)^ 1 and AA = Xm/Vm [2^. 
In addition to these settings dictated by regularization 
issues, we work near the limit when all the components 
of A are equal and continue to large imaginary values, 
i.e., we set |2j| 



A(") =iAo + e 



(a) 



(2.3) 



with real Aq and f satisfying |Ao| 1 and J2a=i '^'"^ ~ 
0. This specific choice for the replica current is tailored 
so that: First, we can assign a sign (positive or negative 
direction) to the multidimensional replica currents. Sec- 
ond, it allows us in actual calculations involving summa- 
tions or integrations over the replica currents to resort to 
the saddle point approximation which is crucial since the 
elementary circuit elements as modelled by Eq. (|2.1|l are 
non-linear. Finally, we invoke the conditions Ag ^ 
and DXl < ^ < 1 making A^ = ^ - DAg a small 
positive quantity so that expansions in powers of A^ are 
reliable. 

After this excursion into formal aspects which we feel, 
however, was necessary for a proper definition of the 



model, we now turn to its physical contents. One of the 
basic quantities describing the electric transport proper- 
ties of DP clusters is the average macroscopic resistance 



Mfl,(x,x') = 



(x+(x,xOi?+(x,xO)c 
P(x,x') 



(2.4) 



when an external current / is inserted at a point x and 
extracted at another point x' provided, that the 2 points 
are positively connected. In Eq. (|2.4|l . _R_|_(x, x') is the 
total resistance if / is inserted at x and extracted at point 
x' and x+(x,x') is an indicator function that gives the 
value 1 one if x and x' are positively connected (if / 
can percolated from x to x'), and zero otherwise. (• ■ ■ )c 
denotes the disorder average over all configurations of 
the diluted lattice and P(x, x') — (x+(x, x'))c is the 
connectivity that measures the probability for x and x' 
being positively connected. The model is set up so that 
the 2-point correlation function 



G2(x,x',A) = (s(x,A)s(x',-A) 



(2.5) 



is a generating function for AIfi[x, x'). This can be seen 
by noting that 



G'2(x, x'. A) = ( exp 



A^ 



i?+(x,x') 



P(x,x')<^l-yAfK(x,x') 



(2.6) 



up to an unimportant multiplicative factor that goes to 
1 in the replica limit D 0. Thus, once we have cal- 
culated the 2-point correlation function, as we are able 
by using renormalized field theory, we can extract the 
average resistance by using 



Mfl(x,x') 



d 



a(-A2/2) 



In G2 ( X, x', A 



A=0 



(2.7) 



B. The nonlinear RRDN 

By generalizing the model we can study much more 
than just the average resistance. Suppose that the di- 
rected bonds obey the non-linear Ohm's law |2^ 



fTfc,. (Vb^. )Vb.. , 1 14,. , 



h 



(2.8) 



The dynamic functional of the nonlinear RRDN consti- 
tuted by these elements is of the same form _as_the func- 
tional H2.2|l . however, with the replacement 



WrAr(A) , 



where r — 1/s and 
A. (A) 



D 

E 

a=l 



-zA(")) 



r+l 



(2.9) 



(2.10) 
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The 2-point correlation function in the generahzed model 
has the property that, up to an unimportant constant, 



G2(x,x', A) = ( cxp 



A.(A) 
r + 1 



i?,.+ (x,x') 



P(x,x')|l-^AfH.(x,x') + 



(2.11) 



where i?r,+ (x, x') is the nonlinear resistance between the 
two terminals. Thus, the 2-point correlation function 
is here a generating function for the average nonlinear 
resistance 



Mi?,_(x,x') 



(X+(x,x')i?r,+ (x,x'))c 



P(x,x') 

which can be calculated by using 
d 



Affl^(x,x') = 



a(-A,,(A)/(r + l)) 



In G2 ( x, x', A 



A=0 

(2.12) 



The generalized RRDN has the benefit that it features 
the free parameter r and that this parameter can be used 
to access several physical quantities. For r ^ 1 one re- 
trieves, of course, the linear average resistance Mi?(x, x'). 
For r , one obtains the mass (the average number 

of bonds) Mb of the backbone, 



Mb (x, x') ^ lim M (x, x') , 

r^ — 1 + 



(2.13) 



as one can see straightforwardly by considering the over- 
all dissipated electric power. Following the lines of Blu- 
menfeld and Aharony 25], one can show that 

Mred(x,x') ~ lim Mi?^(x,x') , (2.14) 

r— !-oo 

for the mass of the red bonds and 

M^in(x,x')- lim Affl,^(x,x'), (2.15) 

for the mass of the chemical path. 



C. The noisy RRDN 

To study multifractal aspects of transport on DP clus- 
ters we can generalize the RRDN so that it features static 
noise. Suppose that the directed bonds have a current- 
voltage characteristic as stated in Eq. (|2.1|l but that the 
value of the conductance of occupied bonds under for- 
ward bias fluctuates statically about some average value. 
To model this effect we set (Jb — ^bjb where <jb is a ran- 
dom variable distributed according to some distribution 
function / with mean ^ — a and higher cumulants A^'-^' 
satisfying A^'^ ^ ct'. The condition on the cumulants is 
imposed to suppress unphysical negative conductances. 



The distribution function / might for example be Gaus- 
sian, however, our considerations are not limited to this 
particular choice. 

In order to access the multifractal properties of this 
noisy RRDN we must treat the average (• • • )c over the 
configurations of the diluted lattice and the noise aver- 
age {•••}/ independently. For this purpose we go beyond 
the usual replica trick and use a second replication pa- 
rameter so that the replicated currents A are {D x E) 
tuples and not merely D tuples |2^. Though the replica 
space for this model is somewhat more complicated than 
the one discussed in Sec. Ill Al essentially the same reg- 
ularization issues appear here and we can make much 
the same settings and impose much the same conditions 
as in Sec. Ill Al albeit in a somewhat more general form. 
For details we refer to Refs. 0, 01 where we derived and 
analyzed a field theoretic model for the noisy RRDN. 
The dynamic functional embodying this model near the 
transitions from the nonpercolating to the directed per- 
colating phases has the same form as the functional (|2.2|l . 
However, we have to make the replacement 0,13 

A^A (2.16) 
and there is an additional term in the dynamic functional. 



1=2 



with dangerously irrelevant interactions 



(2.17) 



(2.18) 



with coupling constants vi ^ A^''/'''^'- Here, the O*-'' are 
composite fields 



(2.19) 



featuring the homogeneous polynomials 

E r D ^-,1 



/3=1 



(2.20) 



The two-point correlation function G2(x,x', A) of the 
noisy RRDN contains a lot of information on the trans- 
port properties of DP clusters. In particular it contains 
the so-called noise cumulants 



C«(x,x') = 



(x+(x,xO{E+(x,xO'}^^^)^ 
P(x,x') 



(2.21) 



where {i?+(x, x')'}^*^"* stands for the ^th cumulant of the 
resistance i?+(x,x') with respect to /. This can be un- 
derstood by noting that, up to an unimportant multi- 
plicative constant that goes to one in the replica limit 
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G2(x, x', A) 



exp 



E 

.1=1 



= P(x,x') 1 + ^ 



(-1/2)' 



1=1 



(2.22) 



i.e., the two-point correlation function is a generating 
function for the noise cumulants which can be extracted 
by using 



Ci;.)(x,x') = 



d 



d 



lnG2(x,x', a; 



A = 

(2.23) 



Although the noise cumulants are interesting in their own 
right, we are primarily interested in a family of observ- 
ables that is more intuitive, viz. the multifractal moments 



M)'^(x,x') 



X+(x,x')E.(V^)' 



c 



P(x,x') 



(2.24) 



of the current distribution. Using Cohn's theorem |27| . 
one can show that 



Mf (x,x')^C)i^(x,x'). 



(0/ 



(2.25) 



This relation allows us to study the multifractal current 
distribution indirectly via the noise cumulants, which is 
important, since, to our knowledge, there is no direct way 
of calculating the moments of the current distribution by 
field theoretic means. 



As usual, a central building block of our renormaliza- 
tion group analysis is a diagrammatic perturbation the- 
ory. Here, the diagrammatic elements for constructing 
Feynman diagrams are the vertex pg and the propagator 



G(p,t,A) = G+(p,i,A) + G_(p,t,A), 



(3.1) 



with 



G±ip,t,X) = ei±t)9{±Xo) 



X exp 



TP 



(r + p2 + WrAriX)) t\ (l - %_o) , (3.2) 



where the factor (1 ^ g) enforces the constraint X ^ 0. 
For the actual calculations it is sufficient to keep ei- 
ther G+(p,i, A) or G_(p,i, A) and we choose to keep 
G+(p,i,A). 

Calculating the conducting Feynman diagrams in di- 
mensional regularization, one encounters ultraviolet di- 
vergences in the form of poles in the deviation s = A — d± 
from the upper critical transversal dimension d± = 4. 
These s poles can be handled by using the renormaliza- 
tion scheme 



s - 

T 

9' 



P - 
Wr 



S = Z ^Zpp, 
Wr = Zp ^Z^ 



= Z-^Z;^Z^G-^up' 



(3.3a) 
(3.3b) 
(3.3c) 



where /i ^ is the usual arbitrary mesoscopic length scale. 
The factor = {4:Tr)-'^^^^r{l + e/2), with F denoting 
the Gamma function, is introduced for convenience. Z, 
Z-r, Zp, and Z^ are the usual DP Z factors known to 
second order in e J29l|30j. In our work on the linear 
RRDN, see Refs. [3, we determined Z^ = Z,,,, to 
second order in e. Studying the nonlinear RRDN 
we showed to two-loop order that 



III. RENORMALIZATION 

In recent years we have investigated the scaling behav- 
ior of the average resistance, the fractal masses, and the 
multifractal moments of the current distribution below 
the upper critical dimension = 4 |^i5|,|^i2i,l3. In par- 
ticular we calculated various scaling exponents and frac- 
tal dimensions by using renormalized field theory. The 
aim of this section is to briefly review central elements 
of our field theory to provide background, to establish 
notation and to gather previous results that we will need 
as input as we go on. For background on renormalized 
field theory in general we refer to p^ . 



Zn = lim Zy_ 



Zr^, 



lim Z^ 
lim Z,„ 



Zr. 



(3.4a) 
(3.4b) 
(3.4c) 



It is not unreasonable to expect that these special rela- 
tions hold to arbitrary order in perturbation theory, i.e., 
that they are Ward identities resulting from the symme- 
tries of the system. Identifying these symmetries and 
proving the presumable Ward identities is a challenging 
open issue for future work. Here, for our calculations to 
follow below we need to know the renormalization factors 
explicitly only to one-loop order. 



A. The linear and the nonlinear RRDN 

To keep our presentation compact, we will review here 
only the general case, i.e., the nonlinear RRDN. The lin- 
ear RRDN can be retrieved by simply letting r 1. 



Z= 1 + — 

4e 



Z,, 



Zt = 1 + — — , Zu = 1 H . 
2e £ 



Z,i,., — 1 



-f 

2£ V 



(3.5a) 
(3.5b) 
(3.5c) 
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The fact that the unrenormahzed theory must be in- 
dependent of the inverse length scale introduced in the 
renormahzation process can be used in a routine fashion 
to set up a Gell-Mann-Low renormahzation group equa- 
tion (RGE) for the correlation functions. For an N point 
function this RGE reads 



V 



N 



G 



(A)};r, 



(3.6) 



where 



V,. 



-a 



du 



TK— + WrCr^ h pCp 

OWr 



dp 



(3.7) 



The functions featured in the RGE are given to two- 
loop order by 



Piu) 

k{u) 
7(u) 
Cp(") 



3m^ / 4\ , 



u 
4 
u 



6 - 9 In ■ 



3J 32 



O (u^) 



Cr(u) =Cr.lU + Cr,2 + O (u^) . 



(3.8a) 
(3.8b) 

(3.8c) 

(3.8d) 
(3.8e) 



The coefHcient of the first order term in Eq. I|3.8a|l is 
known for arbitrary r, 



2r+2 



The coefficient Cr,2 is known only for particular values of 
r. Our two-loop analysis of the linear RRDN with r = 1 
gave 



Ci 



5 

32 



(3.9) 



For the remaining values of r of interest here, (^r,2 can be 
inferred readily from the relations 

C-i{u) ^ j{u) ^ Cp{u) , Coin) = -Cpiu) , (3.10a) 
Coo(w) = «(u). (3.10b) 

stemming from the special relations between the Z- 
f actors. 

In the following we will use for the RGE-functions an 
abbreviated notation of the type f{u) = fiu + f2'i 



For example, we will write Eq. (|3.8a|l as (3{u) — f3iu + 
(32V?' + p^u^ + 0{u^) and likewise for the other RGE- 
functions. 



B. The noisy RRDN 

The parameters vi featured in the dynamic functional 
of the noisy RRDN are, in contrast to the relevant pa- 
rameter w = wi^ dangerously irrelevant, i.e., they are ir- 
relevant on dimensional grounds but they must not be ne- 
glected in studying the noise cumulants because we oth- 
erwise inevitably loose the information we are interested 
in. Due to their irrelevance, the vi cannot be treated 
in the same fashion as the relevant w. A proper treat- 
ment of the vi can be achieved by looking at insertions 
of the dangerously irrelevant interactions V''' into Feyn- 
man diagrams. Due to the irrelevance, insertions of V'''' 
generate a multitude of terms corresponding to interac- 
tions with equal or lower naive dimension than V^'-*. All 
these interactions have to be taken into account in the 
renormahzation process. The interactions of lower naive 
dimension, however, merely lead to subdominant correc- 
tions and can be ignored for our purposes. Keeping all 
the interactions of the same naive dimension, we have a 
renormahzation in matrix form 



(0 



(3.11) 



where V^'-* — (V^'\ Vj''', • • • ) is a vector that contains all 
the interactions generated by V^'^ by the renormahzation 
process including V^'^ itself. The V^'^ are distinguished by 
the feature that the interactions generated by V^'^ do not 
in turn generate , or as we say, that the corresponding 
composite fields O^'^ are master operators j7i. i&. .SL ,32i|. 
Because of their subordinate role, we refer to the remain- 
ing operators Oa'"* and so on as servants. The master 
operators are associated with renormahzation matrices 
Z*-'-* = 1 -|- 0{u), where 1_ stands for the unit matrix, of a 
particularly simple structure, 



^(0 







(3.12) 



\ 

where the * symbolize arbitrary elements. As a conse- 
quence of the simple structure of Z*-'-* the servants can be 
neglected in calculating the scaling index of their master. 
This can be seen as follows. Due to Eq. H2.23|l we are ul- 
timately interested in derivatives of the Green's function 
with insertions of V^'\ G2-\!(i}, with respect to Ki{\) 

evaluated at A = 0. From the renormahzation of the 
master interaction, 



(3.13) 



where the Ya are elements of Z*-'-*, it follows that 



G2(r,i)v(o -Z«G2(r,i) 



v(o+E^a'^G2(r,%(o 

a>2 



(3.14) 
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The coefficients Ya^ pertaining to the servants are re- 
quired to make G2(r, t)v(i) free of e-poles. However, only 
the first term on the right hand side of Eq. H3.14|l gives a 
nonzero contribution when we differentiate with respect 
to Ki{\) and then set A = since only the insertion of 

V*^'^ produces the polynomial structure of Ki{\). Hence, 
as long as we restrict ourselves to the properties of the 
noise cumulants cjj''' (x, x') we only need for our prac- 
tical purposes the multiplicative renormalizations of the 
V*^'\ i.e., we only need the element Z*^'^. We can set all 
the other renormalizations Ya^ pertaining to the servants 
formally to zero. 

Understanding that the renormalization of the danger- 
ously irrelevant interactions V^'-' is subtle but that we can 
ignore these subtleties for our practical purposes, we treat 
their couplings vi in much the same way as w. Doing so 
we have to bear in mind, of course, that this procedure 
only makes sense if we expand all Feynman diagrams in 
powers of vi and truncate this expansion after linear or- 
der (which corresponds to using single insertions) . Then, 
we renormalize the vi by setting 



VI = Z^\i 



(3.15) 



The Gell-Mann-Low function 7*^'^ stemming from Z^'^ is 



7W =7^^ + 7(^)^2 ^0(7.3)^ 



(3.18) 



where 



(0 A 3 



1 



8 22'+i 



(3.19) 



7^') satisfies the relations 7^°-' = C-i = 7 — Cp ^-nd 7^"'^-' 
Ci. 7:^'^ is stated for / = 2, • • • , 5 in Table U 



IV. LOGARITHMIC CORRECTIONS 



Zy^ is known for arbitrary / = 0, 1, 2, • • • to one-loop order 
and for most important values of / to two-loop order. 
Below we will need to know Z„j explicitly to one-loop 
order for calculating the desired logarithmic corrections. 
To this order is related to the Z^^ by 



Zvi — ^W2l-l 



u 
2e 



(3.16) 



The RGE can be solved by the method of character- 
istics whereby one introduces a single flow parameter £ 

TABLE I: The coefficients 7^'^ appearing in Eq. lISTTHll . 





2 


3 


4 


5 


72 


-0.24249 


-0.26523 


-0.26988 


-0.27045 



as we show in the appendix. The two-point correlation 
function is now governed by the RGE 



d_ 

dvi 



xG2(r,pt,wX^;{viKi(x)},T,u,fi^ =0. (3.17) 



and sets up characteristic equations that describe how 
the scaling parameters transform under a change of £. 
The characteristic for the momentum scale fi is particu- 
larly simple and has the solution p,{i) = ^i, i.e., a change 
of £ corresponds to a change of the external inverse length 
scale. With help of the solution to the remaining char- 
acteristics one obtains 



G2 iv,pt,WrAr{X);{viKiiX)},T,u,fi) = {p£Y^Z{£)G2 p£r, {^i£Yp{£)t 



Wr{£)K{X) \vi{£)Ki{\)\ f{£) 



{p£Y 



,u{£),l 

(4.1) 



r 



as a solution to the RGE. This formula applies to the 
nonlinear as well as to the noisy RRDN with the under- 
standing that we have to set = and A = A in the 
former and r = 1 and A = A in the latter case. At this 
stage the scaling solution H4.1|l is still rather formal since 
Z{£), p{£), f{£), Wr{£), v{£) and u{£) require specification. 
Below the upper critical dimension, these quantities dis- 



play power law behavior for ^ — s- described by the well 
known critical exponents of the DP universality class, the 
resistance exponents (f>r, and the multifractal exponents 
ipi of the RRDN. Directly in d± = 4, they depend loga- 
rithmically on £ and hence their behavior is qualitatively 
different from the lower dimensional case. 
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Now we will state and solve the characteristics directly 
for d± — 4. The characteristic for the dimensionless cou- 
pling constant u is given by 



(4.2) 



The solution to this differential equation for e = is 



i = i{u) = iou"'^^/^^ exp 



(4.3) 



where is an integration constant. The remaining char- 
acteristics are all of the same structure, namely 



^d\\iQ{u) 



(4.4) 



where Q is a placeholder for Z, p, r, t/j^, and w/, re- 
spectively, and q is placeholder for 7, Cp, k, Cr, and 7*^'^ 
respectively. Exploiting i d/M — f3d/du one obtains the 
solution 



Q{u) = Qo W 



91 /& 



exp 



(92/^2 - 91/^3) 
PI 



, (4.5) 



where Qo is a non universal integration constant. 

After having reviewed some of the cornerstones of the 
field theory of the RRDN we will now determine the 
critical behavior of the connectivity, the average resis- 
tance, the fractal masses of the backbone, the chemical 
path and the red bonds as well as the multifractal mo- 
ments of the current distribution at the upper critical 
dimension d± = 4. For simplicity, we set in the following 
x' = (r',0) = (0,0) and x = (r,0) = (0,t) and restrict 
our attention to the behavior of the aforementioned quan- 
tities as functions of the time-like variable t. To this end 
we choose 



(4.6) 



where Xq is a constant of order 1. With this choice u 
and £ tend to zero for pn'^t 00. From Eq. (|4.3|l and 



Eq. H4.5|l . specialized to p, we obtain with f32 =3/2 



t = to exp i^—j [1 + 0{u)] (4.7) 

where to = Xo/(/i^^§ po) is yet another non universal 
constant and a is given by 

/32Cp,i-2/33 157 ,53 4 



2/32 



. In^T = 0.97653. (4., 
192 96 3 ^ 



At certain stages it will be more convenient to use, in- 
stead of using the original t, the variable 



s--ln(tAo) 



For s, Eq. 1)4. 7|l translates into 

s = u^^ — alnu + 0{u) 
Finally, we find by using Eq. (|4.10() 



u = s exp 



In s 



O 



' In s In s 1 



(4.9) 



(4.10) 



(4.11) 



for the dimensionless coupling constant as a function of 
the longitudinal coordinate s. 

At this point we know the general scaling form H4.1() 
of the 2-point correlation function and thus we have an 
important part of the information that we need to un- 
derstand the critical behavior of connectivity, the fractal 
masses and the multifractal moments at the upper criti- 
cal dimension. However, for determining the logarithmic 
corrections to the scaling behavior of these quantities be- 
yond the leading correction, we need to know some addi- 
tional information about the scaling function appearing 
on the right hand side of Eq. 14.1|l . 

Expanding the right hand side of Eq. I|4.1|) at criticality 
T = to linear order in Wr^r and viKi, and by using 
Eq. (|4.6|l we find that the two-point correlation function 
is at o?_L = 4 of the form 



G2(0, pt, WrKW]{viKi{\)}, 0, u, fj) = (plfZ{u)G2{Q, Xo, 0; {0}, 0, u, l) 

X |l + {per^Wr{u)Ar{X)g'2{Xo,u) + (^£)-2^«i(M)ifz(A)5^'^(Xo,w) • • • I , (4.12) 

I 



with expansion coefficients g'2 {Xq , u) = 

dlnG2/dwrAr\)^^o and 5^'^(Xo,u) = d\nG2/dviKi\^^f^. 
For the second logarithmic correction we have to cal- 
culate the functions 6*2(0, ATq, 0; {0}, 0, w, 1), g'2{XQ,u) 
and g2\XQ,u) to one- loop order. Some details on these 
calculations can be found in the Appendix. 



A. Connectivity, average resistance and fractal 
masses 

Our one-loop calculation sketched in the Appendix 
yields that the scaling functions relevant for the connec- 
tivity, the average resistance and the fractal masses are 
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given by 

G2 (0, ^0, 0; u, 0, l)=co[l + Apu + 0{u^)] , (4.13a) 
g'^ (Xo, u) ^ci [l + ArU + Oiu^)] , (4.13b) 

where cq and ci are nonuniversal constants and where 
Ap and A^ are amphtudes that can depend on Xq. Ap 
turns out being independent of Xq, 



^--16' 
For the amphtude Ar we find 

1 - In 2 - Z / 1 

= 1 - TT- 

8 V 2"^ 



(4.14) 



(4.15) 



where Z — Z{Xo) is a nonuniversal constant that we 
have defined, to stay consistent with our work on loga- 
rithmic corrections for dynamic properties of DP |14| , as 
Z{Xo) = Ce + ln(2Xo) with Ce = 0.577215... being Eu- 
ler's constant. Alternatively, it is possible to choose the 
arbitrary parameter Xq in such a way that Z — 1 — In 2 so 
that the amplitude Ar becomes simply zero. Physically, 
this would only change the nonuniversal time scale to en- 
tering the scaling functions through Eq. H4.7|l . However, 
we prefer to choose Z as we did for two reasons. First, as 
mentioned above, this allows us to stay consistent with 
previous work. Second, leaving Z in our formulas has the 
advantage that Z can be used as a fit parameter when it 
comes to comparing our results to simulations. 

Now we are finally in the position to assemble our re- 
sults for the connectivity and the average (nonlinear) re- 
sistance. Merging Eqs. H2.11|l . (|4.12|) . and H4.13a|) we 
obtain for the connectivity 



Pit) - i^iif Ziu) [1 + Apu + Oiu^)] 



(4.16) 



Exploiting our choice H4.6|) and using our knowledge 
about the solutions of the characteristics we find 



P{t)/Po 



l + ^)cxp{-cpu) [l + 0{u^)] 



3 \ CP ^ /In^s Ins 1 



(4.17a) 
(4.17b) 



where Pq is a nonuniversal constant and 



cp 



P2 



192 V 



7- 



4 

34 m- 



0.01448. 

(4.18) 



Equation (|4.17a(l in conjunction with Eq. H4.7I) can 
be viewed as a parametric representation of the tuple 
{t^P{t),t) with u serving as the free parameter. Equa- 
tion (|4.17bl) states our result in a more traditional form. 
Though perhaps less intuitive, the parametric form has 
the conceptual advantage that it involves only one ex- 
pansion variable, viz. the effective coupling constant u 



TABLE II: Values of the numbers br and Cr appearing in 
Eqs. K^ . 



r 


-1 


1 


00 


br 


0.16275 


-0.08137 


-0.16275 


Cr 


0.10211 


-0.02519 


-0.03589 



whereas in the traditional form functions of the longitu- 
dinal variable such as 1/s^, Ins/s^, It? s/s^ and so on 
compete against each other. Note that the usually lead- 
ing logarithmic correction is absent in case of the connec- 
tivity P{t). Hence, one has to go to the next order, as we 
did, to see a deviation from the mean-field field behavior 
P{t) cx t-2. 

Next we turn to the average nonlinear resistance. Us- 
ing Eqs. (E21, EH), (|4.13b|) as well as Eq. 
yields 

Mp^it) - {^itj-^ Wr{u)[l + ArU + 0{u'^)\ . (4.19) 

Inserting several intermediate results we are led to 

t-i MR^{t)/MR^^o = [u-^ + B)"^ cxp i-cr u) 

X [1 + O (u2)] (4.20a) 

br In s + ( 



B] 



O 



In^ ,s In s 1 



with a nonuniversal constant Mp^^o Etnd 

Cp.l + Cr.l 1 — 2 ^ 

br = a ar , 

_ (Cp,l + Cr,l)P3 Cp,2 + Cr,2 



Pi 



P2 



B=^ = -^(l-\n2-Z) 
4^ ' 



(4.20b) 

(4.21a) 
(4.21b) 
(4.21c) 

(4.21d) 



As explained above, the physically most important lim- 
its of r are r 1,— l"'",oo,0^ since these limits yield 
respectively the average linear resistance and the fractal 
masses of the backbone, the red bonds, and the chemical 
path. The numerical values of br and Cr in the limits 
r 1,-1 +,oo are collected in Table IHl In the limit 
r ^ 0^ the quantities a^, br, Cr, and Ar vanish and the 
amplitude B does not make much sense. This vanishing 
of the parameters had to be expected because in DP the 
length of the chemical path between the two points (0, 0) 
and {0,t) is essentially t and that hence 



Mr, 



CX t 



(4.22) 



with no logarithmic correction. The fact that our result is 
in conformity with this anticipated behavior is reassuring 
and we rate it as an important consistency check for our 
calculation. 
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TABLE III: Values of the numbers fe^'^ and c'-'\ appearing m 
Eqs. KTEt . 



I 


2 


3 


4 


5 




-0.14241 


-0.15766 


-0.16148 


-0.16243 


c(0 


-0.03263 


-0.03371 


-0.03466 


-0.03530 



B. Multifractal moments 

The logarithmic corrections for the multifractal mo- 
ments can be calculated in much the same way as the 
corrections for the average nonlinear resistance. Our one- 
loop calculation sketched in the Appendix yields 



gf\Xo, u) = ci [1 + A21-1 u + 0{v?)\ , 



(4.23) 



where the amplitude A21-1 can be read off from 
Eg. (jlT^ se tting r = 2 1-1 Using Eqs_g!23> ESS, 
KV^ . as well as ifT^ and Eq. we obtain 

Mf (t) ^ {u) (1 + A21-1 u + 0(u2)) . (4.24) 

Recalling our choice for the flow parameter, Eq. H4.6|l . 
and using that the solution of the characteristic for w^^'^ 
is of the form H4.5II we obtain 



t-^ Mf\t)/M\l = (s-i + Bf" exp(-c(') u)[l + (s^)] 



r(0 



B) 



6(')ln; 



40 



o 



(4.25a) 
In^ s In s 1 

(4.25b) 



where Afj'g ^-"^^ nonuniversal constants. The amplitude 
B is the same as the one in Eq. (14.21(1 . The parameters 
a^') and 6^'' are related to the parameters of Eqs. (|4.21|) 

by 



.(0 



5(0 



02; ^1 



and c^'-* is given by 

(0 _ (Cp,i+7f^)/33 Cp.2 + 



(0 

72 



Pi 



P2 



(4.26) 



(4.27) 



Numerical values for the coefhcient 



,(0 



Table |ini Note that M) 
Note also that Mj' 



(1) 



C 



(0) 



C 



(0) 



since J2bi^b/I) coincides for I 
current-carrying bonds and thus Mj"' = Mb 



can be found in 

= Mr as it should. 

Mb- This must hold 
> with the number of 

(0) 



V. CONCLUDING REMARKS 

In summary, we have studied the electrical transport 
properties of DP clusters by using the methods of renor- 
malized field theory. To our knowledge, logarithmic cor- 
rections for the connectivity, the average resistance, the 



fractal masses of the backbone, the red bonds, the chem- 
ical path and the multifractal moments of the current 
distribution have not been considered hitherto. We cal- 
culated these corrections up to and including the next to 
leading correction. 

We hope that our work triggers complementary numer- 
ical simulations. With todays computer hardware and 
sophisticated algorithms, our results should be testable 
by numerical work. Because we went beyond just calcu- 
lating the leading corrections, we are optimistic that our 
results to compare well with simulations, perhaps even 
quantitatively. 



APPENDIX: AMPLITUDES 

In this Appendix we deliver some details on our calcu- 
lation of the amplitudes entering the second logarithmic 
corrections. First, let us look at the two-point correlation 
function at zero-loop level. As a function of the longitu- 
dinal and the transversal coordinates we have 



(r,t,A) = / exp(zp.r)G+(p,i,A). 



(A.l) 



where is an abbreviation for l/(27r)'*^ J d'^'^p. To treat 
the nonlinear and the noisy RRDN in one go, we use 



G+(p, t, A) = (1 - (5a,o) cxp 



p(T + p- 



L{X))t 



(A.2) 



as our Gaussian propagator. Here, L(A) stands for the 
polynomial 



L(A) = u.,A,(A)-l-^z;/if/(A). 



(A.3) 



As in Sec. IIVI it is understood that we have to set f ; = 
and A = A or r = 1 and A = A to retrieve the nonlinear 
or the noisy RRDN, respectively. 

To one-loop order the correlation function is governed 
by the Dyson-equation 

G2(r,i,A)-Gf (r,i,A) 

exp(zp • r) / dti I dt2G(p,t — ti, X) 
Jo Jo 

X S(p,ti-t2,A)G(p,i2,A), (A.4) 
where E stands for the self-energy S depicted in Fig. [SJ 



S(p,s,A) = 
-h ^exp 



J2 I |-2exp(-pL(A)s) 



p{L{\/2 + k) + L{\/2 - k))i 



y exp [ - p{2t + (p/2 + q)2 + (p/2 - q)')s 



(A.5) 
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FIG. 3: Diagrammatic representation of the self-energy 
E(q, t,A) at one-loop order. 



Note that we have included here an extra term indepen- 
dent of K to switch from the initially restricted current 
summation to a current summation without the restric- 
tion K 7^ 0. The integration over q is straigthforward. 
After switching from the current summation to an inte- 
gration over K we employ the saddle point method which 
we can since Aq 3> 0. From the symmetry of the in- 
tegrand it follows irrespectively of the detailed form of 
L(A) that the locus of the saddle-point is at k = 0. Con- 
sequently, in the limit D — *■ 0, 



S(p,s,A) 



2 



exp(-2/9L(A/2)s) - 2 cxp(-pi:(A)i 



(87rps) exp(-p(2r + pV2)s). 



(A.6) 



One of the time-integrations and the p-integration in 
Eq. I|A.4|1 are easily done. Being interested in critical- 
ity we set r = 0. For simplicity, we restrict our attention 
to r = 0. We expand the right hand side of Eq. (|A.6|I to 
first order in L and u. To prepare for the renormaliza- 
tion step to follow, we mark unrenormalized quantities 
with an open circle, s s, p p and so on as we do in 
the renormalization schemes (|3.3|l and (|3.15l) . Express- 
ing various quantities through their renormalizcd coun- 
terparts as specified by Eqs. I|3.3|l and (j3.15|l we arrive 
at 



2f\2-d±/2 



8r(l + e/2) 



Equation (|3.3|l implies that the correlation function is 
renormalizcd by 

G2(0,t,A) = Z-iG2(0,i,A). (A.9) 

Using again the renormalization factors stated in 
Eqs. (E!ll) and jTTH|l we get 

pL{X) = pWrAriX) + ^ pviKi(X) 
I 

= p(^Z-^Z^,^WrAr{X) + ^ Z-^Z,,viKi{\y^ 
= p(l(A) + J[L(A)-2i(A/2)]), (A.IO) 

where we have used the homogeneity of the polynomi- 
als Ar(A) and Ki{X). From this relation we obtain after 
careful expansion in u and L{X) to first order 



Z-^Gf\o,t,X) 



(1 - m/16) 



d±/2 Y 



(inpt) 

L(A) + J[L(A)-2i(A/2)]]} 



(A.ll) 



Using this result in Eq. (|A.8|I we obtain after a final e- 
expansion 



+ pi-(ln2-l)[L(A)-2i(A/2)] 
■ {2pfiHY^^ 



u 
4e 



[L(A)-2i(A/2)]} 



r(l + e/2 

(1 + 3u/16) 
' {Anpty^ 

pt^ (In 2 - l-f lni2pfih) + Ce) [L{X) - 2L(A/2]) } . 

(A.12) 



{l-p<L(A) 



X / dx- 



[a;(l-a;/2)] 



^jy^ |l + xpt [2L(A/2) - L{X)] I . Taking into account Eq. g^l finally yields 



(A.7) 



The x-integral is performed using dimensional regular- 
ization. We obtain, up to terms of 0(e), 

G.(0,t,A)/Gf(0,,A) = l + gf^ 

2 + pi(2/£ + ln2- 1)[L(A) -2i(A/2)]j . (A.8) 



G2(0,to,A) =co(l + 3ii/16){l-ci 

X {l+Aru)wrAr(X)+Y,{^+^2l-lu)viKi{X) } 

(A.13) 

with nonuniversal constants co and ci and the amplitude 
Ar as stated in Eq. (|4.15|l . 
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